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Based on a simple approximation scheme we have computed the local density of states (LDOS) 
of the antiferromagnetic and ferromagnetic Kondo models for the full range of band occupations 
and coupling strengths. For both models the LDOS with its full energy dependence has not been 
calculated before. Arguments are given for the results to be qualitatively trustworthy despite the 
simplicity of the approximation scheme. 



I. INTRODUCTION 

The term antiferromagnetic Kondo model stands for 
the well-known single-impurity Kondo model, whose ori- 
gin is ascribed to ZcncrEJ. The attribute antiferromag- 
netic is used to stress the difference to the ferromagnetic 
Kondo model. In both models the spin of a magnetic im- 
purity is coupled to the spins of the conduction electrons 
of a non-magnetic host lattice. In the (antiferromag- 
netic) Kondo model the sign of the coupling constant is 
such that antiparallel coupling between the impurity and 
the conduction-electron spins is favoured (antiferromag- 
netic intraatomic exchange) . In the ferromagnetic Kondo 
model, with the opposite sign of the coupling constant, 
parallel coupling is favoured (ferromagnetic intraatomic 
exchange) . 

The (antiferromagnetic) Kondo model has been one of 
the most intensively discussed many-body models in solid 
state physics since Kondo in 1964 succeeded in explaining 
the resistance minimum of metals with small amounts of 
transition element impurities 2 !. Thorough theoretical in- 
vestigations later on led to the discovery of the cause of 
the minimum, the Kondo effect, which is the collective 
(many-body) screening of the impurity spin by the spins 
of the conduction electrons below the Kondo temperature 
Tr-lIu- By means oOVilson's numerical renormali^ation 
group theory (NRG)Eil and Bethe-ansatz methodstJu exact 
solutions for thermodynamic properties such as the mag- 
netic susceptibility and the heat capacity were obtained. 
Later the NRG was successfully extended to the^mou- 
tation of dynamic (energy-dependent) quantitiea3tl'E£l. 

Although the antiferromagnetic Kondo model can thus 
be regarded as in principle solved, its local density of 
states (LDOS) with the_full energy dependence has not 
been calculated so faiO. Nevertheless the LDOS is 
an interesting quantity, particularly since scanning tun- 
neling spectroscopy (STS) has made it experimentally 
accessibleE3lSolia. 

The LDOS of the ferromagnetic Kondo model, too, 
has not been calculated yet. The ferromagnetic Kondo 
model has generally received less attention than its anti- 
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ferromagnetic counterpart. This is probably due to the 
fact that a spectacular "ferromagnetic Kondo effect" does 
not exist because the effective coupling does not scale to 
infinity as in the antiferromagnetic Kondo modelEKa. A 
second reason might be that there are fewer experimental 
realizations than in the antiferromagnetic case. 

In this paper we present results for the T = 0-LDOS of 
both the antiferromagnetic and the ferromagnetic Kondo 
models for the full range of band occupations and cou- 
pling strengths, gained by means of a simple Green's 
function approximation scheme based on Nagaoka's de- 
coupling procedureEa. The equations of Nagaoka's de- 
coupling scheme were solved analytically by Zittartz and 
Mullcr-HartmanntZl. However, their equations still de- 
pend self-consistently on the scattering matrix. There- 
fore, in order to find self-consistent numerical solutions, 
it is equally well justified to use Nagaoka's original equa- 
tions. 

Qualitatively results prove trustworthy despite the 
crudness of the approximations. For weak couplings the 
main correlation feature in the impurity-site LDOS of 
the antiferromagnetic Kondo model is an antiresonance, 
which is slightly shifted with respect to the Fermi level. 
There is an analogous antiresonance in the LDOS of 
the single-impurity Anderson model. A related antires- 
onance structure was observed in re cent i ST S measure- 
ments on magnetic impurity adatomsr 3 IF 4 rl 15 l. For strong 
couplings the dominant characteristic are quasiparticle 
(5-peaks. These can be traced back to the limiting case of 
an empty conduction band, for which exact results and 
rigorous interpretation are given, and the atomic limited. 

In the ferromagnetic Kondo model the LDOS at the 
Fermi level is practically independent on the coupling 
strength, which is consistent with the fact of a vanish- 
ing phase shift. However, the LDOS is diminished in the 
vicinity, giving rise to a peak-like structure. If represent- 
ing a true feature, this peak structure should be mea- 
surable in STS experiments on magnetic adatoms with 
ferromagnetic intraatomic exchange. 

It is clear that the approximation scheme is too sim- 
ple to yield quantitatively reliable results. Subtle features 
like the exponential Kondo scale in the antiferromagnetic 
Kondo model cannot be reproduced. Much better meth- 
ods like numerical renormalization|-eroup theory (NRG) 
or Quantum Monte Carlo (QMC)ll3 would be required 
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for that. 

This paper is structured as follows: First the exact so- 
lution for the limiting case of an empty conduction band 
(Sec. HA) a nd t he approximation scheme for finite band 



occupation ( II B ) are explained. Then the results for the 
LDOS are discussed in detail: first, for the case of an 



empty conduction band ( [II A ) , which is similar in both 
models, second, for the case of finite band occupation 
in the antiferromagnetic Kondo model ( III B ), and third, 
for finit e ban d occupations in the ferromagnetic Kondo 
model ( III C ) . After that the influence of the magnetic 
impurity on occupation numbers is discussed with em- 
phasis on a seemingly contradictory result by other au- 
thors ( III D ). Finally we draw some conclusions (IV). 



The free Green's functions are given by g[^\e) 
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Solving (H) - (||) yields the solution: 
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With (|g), the LDOS pj\ of any given site Rj can be cal- 
culated: Pj]{E) = — ilmG^i (E). The spin dependence 
of Pj-f is purely formal. 

For our calculations we have chosen the tight-binding 
dispersion with nearest-neighbour hopping for the s.c. 
lattice: 



II. THEORY 

The Hamiltonian of the antiferromagnetic and ferro- 
magnetic Kondo models is given by 
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Replacing <x by Fermi operators yields 



e(k) = — 2t(cos k x a + cos k y a + cos k z 



(7) 



where t is the hopping integral and a the lattice con- 
stant. (Q) implies a symmetric Bloch density of states 
(BDOS) and ImG^-E) = ±ImG ( i ) (~E), which result 
in a mirror symmetry of the LDOS when changing from 
the antiferromagnetic to the ferromagnetic Kondo model: 
pffi(E) = p^ J \—E) at any site R,. 



B. Approximation scheme for finite band 
occupation 



with z a = 5 a] - 5 al and Sg = 5 a1 S£ + S al S Q . The 
first part of H with the hopping integrals describes 
the free movement of conduction electrons in a nonde- 
generate band. c^J annihilates (creates) an electron in 
a spin-cr Wannier state at the lattice site R^. The sec- 
ond part stands for the interaction between the impurity 
spin So and the spins of the conduction electrons repre- 
sented by <r . Throughout this paper a spin-i impurity 
will be considered. We assume a ^-independent coupling 
constant J. Therefore the coupling is effectively local, 
involving only electrons at the impurity site Ro- TL rep- 
resents the (antiferromagnetic) Kondo model if J < 0, 
and the ferromagnetic Kondo model for J > 0. 



A. Exact solution for an empty conduction band 

Appropriate Fourier transformations of the equa- 
tions of motion for the Green's functions G^^E) = 

(( c iT ; c] T )) E and T 0in (E) = (( Sfic^ ; c] T )) £ lead to the 
following closed system of equations: 
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For G^ Eqs. and (||) remain correct. For T, 
new Green's functions come ijatp play, to which Nagaoka's 
decoupling scheme is appliedll3. Consequently, Eq. (|^) is 
replaced by 
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We have defined special energy-dependent annihilation 
operators 
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ViV ^ E - e(k) + i0- 
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to avoid multiple fc-space summations within the self- 
consistency cycle of this approximation scheme. The re- 
maining problem is to determine the correlation functions 
of Eq. (|^). They are complex quantities, whose real and 
imaginary parts must be determined separately via the 
spectral theorem with the help of appropriately defined 
operators and Green's functions. For example: 



Re(c 0T c 0T (£;)) = (c 0T c QT (£)) 
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V meaning "principal part of" . One needs 

G 00 ^E,E')=((c m (E);c m )) E > (12) 

with E a parameter and E' the usual energy variable of 

the Green's function to determine (cq^c q1 -(E)) via the 
spectral theorem 
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Comparing the respective equations of motion for 

n r 

G 00 ^{E, E') and G 00 ^(E) one can show that Gqo-\(E, E) 
is fully determined by the "simple" one-electron Green's 
function through 
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Thus Re( c 01 -c Q ^(E) ) and, analogously, Im( cLc Q |(-B) ) 
are fully determined by G 0Q +. 

To determine the real part of the other correla- 
tion function, Re( SgcLe 01 .(-E) ), one needs the Green's 



function r 000T {E, E') = (( c 0T (£) ; c 0T )) B > ■ We as- 
sume the same proportionality (Q) that holds between 

R 

G 00 ^{E, E') and G Q0 ^(E') for the T-Green's functions: 
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Note that Eq. ( |l6| ) is required to work only within 
the spectral-theorem integration over E' to determine 
c ot^ot(^'))- Within that integration great contri- 
butions arise for E' w £7, for which ( |l6| ) can be shown to 
be reasonably justified. Analogous considerations can be 
made for Im( SqCq+Cq*(E) ). 

Hence, one has a closed system of equations consis ting 
of ®, d§, d), ©, ©, the analogues to (§1) and @ 
for determining Im( Cg^c ^.(-B) ) and Im( 5qCq^Cq|(-E) ), 
and the spectral theorems (such as (|l3|)) for the real and 
imaginary parts of the correlation functions. The LDOS 
Pif(E) = —^ImG u ^(E) can be self-consistently deter- 
mined for any lattice site R^. The spin dependence of 
Pjl is purely formal. The occupation number at Rj is 
gained by integrating over the LDOS: 
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III. RESULTS AND DISCUSSION 



A. LDOS in the limiting case of an empty 
conduction band 
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FIG. 1: LDOS of the antiferromagnetic Kondo model at the 
impurity site in the limiting case of an empty conduction band 
for various values of the coupling constant J. 5-peaks are rep- 
resented by pins whose heights are proportional to the corre- 
sponding spectral weights. The free Bloch dispersion is given 
by Eq. (M). The energy is measured in units of the bandwidth. 

We first discuss the limiting case of an empty conduc- 
tion band, in which results are exact and allow rigorous 
interpretation, and which is similar in the antiferromag- 
netic and ferromagnetic Kondo models. The LDOS rep- 
resents the outcome of an inverse-photoemission exper- 
iment in which one test-electron is put into the empty 
conduction band. It thus provides information on the 
impurity's effects on one-electron states. In Figs, [j] and 
^ the LDOS of the antiferromagnetic Kondo model at 
the impurity and the nearest-neighbour sites is shown 
for various coupling constants J < 0. The LDOS of 
the ferromagnetic Kondo model (J > 0) is simply given 
by the mirror image of the LDOS of the antiferromag- 
netic Kondo model with the same |J|. This symmetry 
is a consequence of the Bloch dispersion we have cho- 
sen (Eq. (0)). A general dispersion gives similar but not 
exactly symmetric results for the two models. 

At J = the LDOS corresponds to the free Bloch 
density of states (BDOS). Small couplings cause defor- 
mations of the LDOS, most pronounced at the impu- 
rity site. For J < —0.44 there appears a S-peak below, 
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FIG. 2: LDOS of the antiferromagnetic Kondo model at a 
nearest-neighbour site of the impurity in the limiting case of 
an empty conduction band. The rest as in Fig. |l|. 



for J < —1.33 a second <5-peak above the quasiparticle 
band. The upper peak has greater spectral weight than 
the lower. At the nearest-neighbour site the weight of 
both peaks is considerably less than at the impurity site. 
With increasing \J\ the peaks move away from the quasi- 
particle band. At the impurity site their weight increases, 
while the weight of the quasiparticle band gradually van- 
ishes. At the neighbouring site the 5-peaks lose weight 
(not recognizable in Fig. |2|) and the quasiparticle band 
approaches a certain asymptotic form. 

This can be well understood from the strong-coupling 
limit ( J — > — oo). In that limit the singlet and triplet 
states with the electron completely localized at the im- 
purity site are energy eigenstates. The other one-electron 
states are extended band states, deformed in the vicin- 
ity of and zero at the impurity site, thereby meeting the 
condition of orthogonality to the fully localized states. 
Therefore, in the strong-coupling limit the LDOS at the 
impurity site consists of two 5-peaks, their positions tend- 
ing towards — oo and +oo, corresponding to the diverg- 
ing energies of the bound singlet and triplet states. The 
quasiparticle band has zero weight as there is no overlap 
of the band states to the impurity site. At the other sites 
the LDOS vice versa consists of a quasiparticle band only. 
It is deformed because the probability of the electron to 
stay at the respective site is redistributed between and 
within the extended states. 

For finite J the bound states have finite overlaps to 
the neighbouring sites of the impurity, which is reflected 
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FIG. 3: Quasi-exponential decrease of the spectral weights 
of the singlet and triplet peaks with increasing distance from 
the impurity site (0,0,0). The weight of the singlet peak is 
shown for \J\ — 1, the weight of the triplet peak for \J\ =2.2. 
(Calculated for same Bloch dispersion as in Fig. 



by quasiparticle ^-peaks with small but finite weights. 
Since the eigenstates with one electron can be classified 
as pure singlet or triplet states for all J, the lower and 
upper quasiparticle <5-peaks correspond to bound singlet 
and triplet states, respectively. The greater weight of 
the upper peak ("triplet peak") results from the degen- 
eracy of the triplet states. The degree of localization of 
the bound states around the impurity site is very large 
as can be seen when comparing the spectral weights of 
the singlet and triplet peaks at the impurity and neigh- 
bouring sites (Fig. |3|). There is a quasi-exponential de- 
crease of weights/overlaps with increasing distance from 
the impurity. The extended band states in case of finite 
J have a finite overlap to the impurity site, which is re- 
flected by the finite weight of the quasiparticle band in 
the impurity-site LDOS. 

If considering singlet and triplet states separately, the 
situation is similar to the quantum-mechanical problem 
of a particle in a potential landscape with a well or bar- 
rier. For instance, the bound singlet state corresponds 
to a bound state in a potential well, which has an ex- 
ponentially decreasing overlap into the adjacent region if 
the depth of the well is finite. In case of an infinitely 
deep well, corresponding to the strong-coupling limit, 
each bound state is restricted to the well while the ex- 
tended states have no overlap into the well region. 

For large values of J the effect of the hopping from 
and to the impurity site becomes negligible compared to 
the spin-spin coupling. Therefore, as far as the impurity 
site is concerned, a comparison with the atomic limitEJ 
is sensible. In case of zero band occupation the LDOS 
in the atomic limit consists of two quasiparticle 5-peaks 
whose positions are at +|j (singlet) and — |J (triplet), 
having spectral weights of j and | , respectively. In Fig. |] 
it is shown how the weights and positions of the singlet 
and triplet peaks in the impurity-site LDOS converge for 
large J towards the corresponding atomic- limit values, 
which confirms the emergence of fully localized (atomic) 
states in the strong-coupling limit. 
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FIG. 4: Convergence of the weights and positions of the 
impurity-site singlet and triplet peaks toward the respective 
atomic-limit values. The solid lines in the insets show the rel- 
ative positions of the peaks (ratio of position in bandwidth=l- 
case and atomic limit). (Calculated for same Bloch dispersion 
as in Fig. El.) 



B. LDOS in the antiferromagnetic Kondo model 

In Fig. ||our results for the Kondo-model LDOS at the 
impurity site for different band occupations n and the 
full range of weak to strong couplings J are presented. 
At more than half filling (n > 1) the LDOS is simply 
given by the mirror image of the corresponding LDOS 
at a band occupation of 2 — n (due to particle-hole sym- 
metry of the Kondo Hamiltonian (Q) and the selection of 
a symmetric free BDOS). For weak to intermediate cou- 
plings (|J| < 0.4) and not too small band occupations 
(n > 0.2) the most prominent feature in the LDOS is 
an antiresonance at the Fermi level, which will be dis- 
cussed in subsection 1. In the strong-coupling regime 
(| J | > 0.6) quasiparticlc <5-peaks are the dominant char- 
acteristic, which are discussed in subsection 2. 



1. Kondo antiresonance 

For small couplings (J = — O.f) the LDOS exhibits a 
narrow, dip-like structure at the Fermi level. When in- 
creasing the coupling strength (J = —0.2, —0.4) it deep- 
ens and broadens, developing into a marked antireso- 
nance. As shown in Fig. ^j, the antiresonance appears not 
exactly at the Fermi energy but is slightly shifted (up- 
ward shift for n < I). For strong couplings (\J\ > 0.6) 
the LDOS at the Fermi level is practically zero. However, 
no antiresonance in the narrower sense is left, instead the 
LDOS is dominated by quasiparticle peaks. In the fol- 
lowing paragraphs we will discuss independent evidence 



for an antiresonance in the (antiferromagnetic) Kondo 
model. After that, the antiresonance will be related to a 
structure recently observed in STS experiments. 

In Ref. ^0] Mezei and Zawadowski considered the 
Kondo model for different types of scattering. They made 
explicit statements on the LDOS for the resonance en- 
ergy, for which they could refer to a definite value of 
the phase shift. For s-type scattering they found com- 
plete suppression of the LDOS, while non-s-type scatter- 
ing was shown to involve no changes at the resonance 
energy. As each of the scattering types (s,p, etc.) is 
connected with a special wave- vector dependence in the 
coupling constant, it seems a priori unclear how to ap- 
ply those results to the Kondo model with a wave-vector 
independent coupling-constant J discussed in this paper. 

Evidence for an antiresonance in the LDOS of the 
Kondo model is provided by a comparison jssrith re- 
sults for the (single-impurity) Anderson modelEH. As is 
well-known, the Anderson model in the so-called Kondo 
limit can be mapped on the Kondo -model by the uni- 
tary Schrieffer- Wolff transformatiorcj. A correspon- 
dence must be expected between results in the Kondo 
limit of the Anderson model and the weak-coupling 
limit of the Kondo model. In the Anderson model 
one does observe an antiresonance in the conduction- 
electron LDOS. It is closely related to the appearance 
of the well-known Kondo resonance (peak) in the impu- 
rity quasiparticle density of states (impurity QDOS). Sol- 
lie and Schlottmann calculated the conduction-electron 
LDOS and the impurity QDOS of the single-impurity 
Anderson model within the slave-boson approached The 
conduction-electron LDOS shows the antiresonance ap- 
pearing at the same energy and with the same width 
as the Kondo resonance in the impurity QDOS, which 
reveals the intimate relation between both structures. 
Both appear not exactly at the Fermi energy but slightly 
above, similar to what we observe for the antiresonance 
in the Kondo model. More recently Schiller and Hersh- 
field emphasized the direct correspondence between an- 
tiresonance and Kondo resonance-.("Abrikosov-Suhl res- 
onance" ) in the Anderson modclEil. 

We have performed a calculation for the Anderson 
model within the modified perturbation theory (MPT). 
The MPT is an extension ofr-tihe conventional second- 
order perturbational schcmeuuS to improve the strong- 
coupling behaviour away from the symmetric point, 
which is the only point at which the conventional pertur- 
bation .theory also converges for lar ge H ubbard interac- 
tion U£3 For details we refer to Ref. £8|, where the MPT 
was successfully applied to the single-impurity Anderson 
model. It has also been used in the context of the dynam- 
ical mean-field theory (DMFT) to examine-thii Hubbard 
modelEj and the periodic Anderson modelEJ'Eil. Our re- 
sults for the single-impurity Anderson model for three 
different band occupations n are shown in Fig. [t]. The 
Kondo resonance in the impurity QDOS and the antires- 
onance in the conduction-electron LDOS are clearly to be 
seen at the Fermi level (for n = 0.2, 0.6 slightly above). 
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FIG. 5: LDOS of the antiferromagnetic Kondo model at the impurity site for various coupling constants J and band occupations 
n. The position of the chemical potential fi is indicated by the thin dotted line. The narrow resonances outside the quasiparticle 
band are quasiparticle <5-peaks. (Calculation for the same Bloch dispersion as in Fig. |l|.) 
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FIG. 6: Antiresonance in the impurity-site LDOS of the anti- 
ferromagnetic Kondo model (n = 0.6, J = —0.1). The Fermi 
level is indicated by the thin dotted line. The shift of the an- 
tiresonance relative to the Fermi level is of the order of 10~ 3 
in units of the bandwidth. (Same calculation as in Fig. |^.) 



With the parameters chosen, an approximate correspon- 
dence to the Kondo model with J = —0.16 is to be ex- 
pected according to the Schrieffer- Wolff transformation. 
For each of the three band occupations the Anderson- 
model LDOS is indeed quite similar to the Kondo-model 
LDOS for J = —0.2 in Fig. |5[ The only basic difference is 
for the symmetric Anderson-model LDOS, where the an- 
tiresonance goes right down to zero, whereas in the n = 1 
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FIG. 7: Conduction-electron LDOS at impurity site (solid 
line) and impurity QDOS (broken line) in the Anderson model 
for different band occupations n calculated with the modified 
perturbation theory (MPT). The position of the Fermi level 
is indicated by the thin dotted line. Parameters: U = 1, 
V — 0.2 (Bloch dispersion as in Fig. o). 



Kondo-model LDOS it does not. For the symmetric An- 
derson model one can show, using the fact that the elec- 
trons form a Fermi liquid, that the conduction-electron 
LDOS of the impurity site must be zero at the Fermi 
energy. Assuming this is conserved in the Schrieffer- 
Wolff transformation, it is apparently a shortcoming of 
our approximation scheme that for small to intermedi- 
ate couplings the symmetric Kondo-model LDOS is not 
fully suppressed at the Fermi energy. Furthermore, one 
cannot expect our approximation scheme to satisfy the 
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subtle exponential dependence of the width of the an- 
tiresonance on the coupling constant J via the Kondo 
temperature Tr-Q, which should occur in analogy to the 
Anderson model. Nevertheless, the antiresonances com- 
ing out from the calculations clearly indicate the correct 
antiresonance of the Kondo model. Since the antircs- 
onance of the Kondo model is closely connected with 
the one in the Anderson model and because of the in- 
timate relation of the latter to the well-known Kondo 
resonance, we use for both antiresonance structures the 
term "Kondo antiresonance" . 

Recently experimental evidence was found for the 
Kondo antiresonance in the conduction-electron LDOS 
of magnetic impurity adatoms on metal surfaces. In 
scanning tunneling spectroscopy (STS) experiments the 
differential conductance dl/dV was measured for the 
cases ofpCje adatoms on a Ag(lllXsurfaceE3, for Co on 
Au(lll)EJ and for Co on Cu(lll)li3. In each experiment 
a narrow antiresonance structure with a slight upward 
shift relative to the Fermi level was observed in the dl/dV 
spectrum. Following Ref. in case of weak tunneling, 
low enough temperatures (T — > 0) and in an energy range 
in which the LDOS of the STM tip is nearly constant, the 
differential conductance dl/dV is in general proportional 
to what we shall call a "mixed LDOS". Contributions 
to this "mixed LDOS" come from all electronic states 
( "tunneling channels" ) of the substrate that are involved 
in the tunneling process. The "mixed LDOS" is not just 
the sum over partial densities of states but also contains 
quantum-interference effects between the different tun- 
neling channels. 

As for their measurements on Ce/Ag(lll), Li et al. as- 
sumed that they wer£-."locally sensitive to the hybridized 
sp conduction band"E3. This is supported by Schiller and 
HershfieldcJ, who achieved optimal theoretical agreement 
with the measured antiresonance curve on the assump- 
tion that alternative tunneling into the impurity (Ce) 
/-orbital was zero. Accordingly, the pure conduction- 
electron LDOS was probed, which allows to directly iden- 
tify the antiresonance in the dl/dV spectrum with the 
Kondo antiresonance in the conduction-electron LDOS. 

By contrast, in the case of Co/Au(lll) the involvement 
of the tunneling channels and the interpretation of the 
dl/dV spectrum are a matter of controversy. Madhavan 
et al. assumed significant contributions from tunneling 
both into the conduction band and the impurity (i-orbital 
of the Co adatom, attributing an asymmetry of the an- 
tiresonance to quantum-interference effects betssieen the 
d-orbital and the conduction electron channelaiJ. This 
was confirmed by an excellent theoretical fit of the an- 
tiresonance. However, Ujsaghy et al. assumed tunneling 
into the conduction-band only and achieved an equally 
good theoreticaLfitnj. The ambiguity is already inherent 
in Fano's theoryLj, which in Refs. 14 and 09 was general- 



ized to the interacting case. For non-interacting electrons 
one can show that the asymmetric shape of a "Fano res- 
onance" depends, first, on the ratio of tunneling into the 
impurity orbital and tunneling into the conduction band, 



and secondly, on the symmetry/asymmetry of the free 
conduction-electron BDOS and the Fermi-level position. 
Each of these factors and combinations of the two can 
lead to the same line shape. This seems to hold in the 
present (interacting) case, too. We conclude that in the 
case of Co/Au(lll) the dl/dV spectrum possibly repre- 
sents a direct measure of the conduction-electron LDOS. 
If so, the asymmetric antiresonance (Fano resonance) di- 
rectly represents a Kondo antiresonance. If not, a Kondo 
antiresonance may still be considered as the main under- 
lying structure apart from impurity-orbital contributions 
and quantum-interference effects. 



2. Quasiparticle 8 -peaks 

In the strong-coupling regime (\J\ > 0.6) the LDOS of 
the Kondo model is dominated by quasiparticle <5-peaks. 
In Fig. ^jthey appear as narrow resonances. For t he limit- 
ing case of an empty conduction band (Sec. Ill A ) we have 
discussed the emergence of a singlet and a triplet peak, 
corresponding to the excitation of bound one-electron 
singlet and triplet eigenstates through inverse photoe- 
mission. With increasing J first the singlet peak emerges 
from the quasiparticle band. In Fig. |B| this can be seen 
again in the first column on the left, which represents the 
limiting case of an empty conduction band. If the singlet 
peak is present, one has to distinguish another case of 
zero band occupation n = (second column in Fig. |^), 
characterized by the chemical potential ji lying between 
the singlet peak and the quasiparticle band. In this case 
there is exactly one electron in the system, occupying 
the bound singlet state, before the excitation in inverse 
or direct photoemission. With this electron the band 
occupation n (i.e. the average number of electrons per 
lattice site) is still zero, but the physical situation differs 
from "the limiting case of an empty conduction band" . 
This is reflected by the drastic differences between the 
corresponding LDOS's in Fig. ^. 

For certain values of the coupling constant (cf. J = 
—0.6) the emergence of a second quasiparticle peak can 
be triggered by increasing the band occupation. The 
generic situation of the strong-coupling regime is the 
presence of two quasiparticle peaks. In the case of half 
filling they are positioned symmetrically above and below 
the quasiparticle band, whose weight gradually vanishes 
with increasing coupling strength (cf. J = —0.6,-0.8). 
The LDOS for large J approaches the LDOS of the 
atomic limit, which for n = 1 consists of two symmet- 
ric singlet peakal3, the lower one observable in direct, 
the upper in inverse photoemission. The quasiparticle 
peaks of the strong-coupling regime clearly correspond 
to the singlet peaks of the atomic limit. Nevertheless, in 
the present case of finite hopping, admixtures of triplet 
coupling at the impurity site are to be expected, which 
only disappear in the limit J — > oo. Note that these 
would not contradict the overall singlet nature of the 
Kondo-model ground state. In Ref. 23 similar quasipar- 
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FIG. 8: LDOS of the ferromagnetic Kondo model at the impurity site for various coupling constants J and band occupations 
n. The position of the chemical potential is indicated by the thin dotted line. The narrow resonances outside the quasiparticle 
band are quasiparticle 5-peaks. (Calculation for the same Bloch dispersion as in Fig. |l|.) 



ticle peaks close to the band edges were observed in the 
conduction-electron LDOS of the single-impurity Ander- 
son model. 



C. LDOS in the ferromagnetic Kondo model 

Figure || shows our results for the impurity-site LDOS 
in the ferromagnetic Kondo model for low to strong 
couplings J and different band occupations n. Due to 
particle-hole symmetry, it is sufficient to consider band 
occupations up to half filling. The most remarkable fea- 
ture is a peak structure at the position of the Fermi level. 
Note that the peak structure is not a true resonance — 
in the sense that the LDOS is not really enhanced at the 
maximum in comparison to the free BDOS. This agrees 
with what one should expect because of the vanishing 
phase shift at the Fermi level in the ferromagnetic Kondo 
model. Instead, compared to the BDOS the LDOS is di- 
minished in the vicinity of the maximum position. The 
peak structure is slightly shifted relative to the Fermi 
level (Fig. ^|) similar to the Kondo antiresonance of the 
antiferromagnetic Kondo model. 

If the peak structure at the Fermi level is a true fea- 
ture of the ferromagnetic Kondo model, it should in prin- 
ciple be detectable in an STS experiment similar to the 



Kondo antiresonance (cf. Sec. IIIB1). In an STS mea- 



surement the differential conductance dl/dV represents 
a direct measure of the conduction-electron LDOS on 
condition that tunneling occurs only into the conduction 
band. The peak structure in the LDOS of the ferromag- 
netic Kondo model should thus be observable as a seem- 
ing enhancement in the differential conductance near the 
Fermi level. However, the preparation of adatoms with 
ferromagnetic exchange (J > 0) is required. Gd adatoms 
might be suitable candidates since in bulk Gd the in- 
traatomic exchange is ferromagnetic, which can be in- 
ferred from an excess T = magnetic moment due to 
induced "ferromagnetic" spin-polarisation of 5c? — 6s con- 
duction electronsEj. 

As in the antiferromagnetic Kondo model, for suffi- 
ciently large couplings J quasiparticle <5-peaks appear, 
whose emergence can partly be triggered by increasing 
the band occupation (cf. J = 0.1). Again, the generic 
situation in the strong-coupling regime is the presence of 
two quasiparticle (5-peaks. In case of half filling they are 
located symmetrically below and above the quasiparticle 
band, whose weight gradually vanishes when increasing 
the coupling strength. Similar to the antiferromagnetic 
Kondo model and the limiting case of an empty con- 
duction band, the LDOS in the strong-coupling regime 
should approach the LDOS of the atomic limit. The two 
symmetric quasiparticle (5-peaks therefore-correspond to 
the two triplet peaks in the atomic limitta of the ferro- 
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FIG. 9: Peak structure in the impurity-site LDOS of the ferro- 
magnetic Kondo model (n = 0.6, J = +0.2). The Fermi level 
is indicated by the thin dotted line. The shift of the peak 
structure relative to the Fermi level is of the order of 10 -3 in 
units of the bandwidth. (Same calculation as in Fig. ^.) 

magnetic Kondo model for n = 1. Nevertheless, in the 
case of finite hopping, admixtures of singlet character are 
to be expected, which only disappear in the limit J — > oo. 
The assumption of singlet admixtures is supported by the 
observation that for J = 1.0 the upper peak in the n = 1- 
LDOS continuously evolves from the pure singlet peak of 
the n = 0-LDOS (limiting case of an empty conduction 
band) when changing the band occupation from to 1 . 

D. Occupation numbers 




-0.5 -0.25 0.25 ,. 0.5 



FIG. 10: Difference between the occupation number at the 
impurity site and the band occupation, (no) — n, in depen- 
dence on the Fermi level jj, for different values of the coupling 
constant J. fj, — —0.5 corresponds to zero band-occupation, 
(i = to half filling and a = +0.5 to a completely filled band. 
(Calculated for same Bloch dispersion as in Fig. Ill) 



In Fig. [L0| the effect of the magnetic impurity on the 
occupation number (no) at the impurity site is demon- 
strated. Depending on the band occupation n, the impu- 
rity can effectively attract or repel the conduction elec- 
trons both in the antiferromagnetic and ferromagnetic 
Kondo models. For a symmetric BDOS the dependence 
on the band occupation is exactly symmetric: below half 
filling the impurity effectively attracts the conduction 
electrons, (no) — n > for —0.5 < fi < 0, above half 
filling it repels them, (n ) — n<0for0<«< 0.5. 

The physical reason for this behaviour is clear and 
should analogously hold for a non-symmetric BDOS. One 
of the possible couplings of an electron to the impurity 
spin (in the antiferromagnetic Kondo model the singlet, 
in the ferromagnetic Kondo model the triplet) decreases 
the energy of the system, making it favourable for the 
electron to stay at the impurity site. Thus, (no) is in- 
creased compared to the band occupation n in the case 
of less than half filling. However, it is unfavourable for 
a second electron to be at the impurity site (since it 
would break the favourable coupling of the other elec- 
tron), which leads to an occupation number less than 
the band occupation in the case of more than half filling. 
In a generalized form this argument implies a critical 
band occupation above which the occupation number at 
the impurity site is less than in the rest of the lattice, 
corresponding to a repulsive effect of the impurity on the 
conduction electrons. Below we discuss results of Refs. |34| 
and which are seemingly contradictory to this. 



The effect of the magnetic impurity on the occupation 
numbers (rij) at the surrounding lattice sites is demon- 
strated in Fig. [fl| Oscillations around the value of the 
band occupation which decrease in amplitude with grow- 
ing distance from the impurity are to be seen. The os- 
cillations in the occupation number reflect the Friedel 
oscillations in the charge density, which are the typical 
reaction of an electron system to an impurity potential. 
Everts and Ganguly derived an expression for the charge 
density in the Kondo model (Ref. |34| ). Sokcevic, Zlatic 
and Horvatic calculated charge-density oscillations in the 
Anderson model (Ref. |3^). A quantitative comparison 
between charge-density and occupation-number results 
is only partly possible. The occupation number is the 
number of electrons in the Wannier states at a given lat- 
tice site, whereas the charge density corresponds to the 
number of electrons per unit volume for a given point of 
space. In contrast to the occupation number, the charge 
density at and near a site depends also on electrons in 
Wannier states of other sites because in real space Wan- 
nier wave-functions are not confined to "their" sites. Al- 
though there is no rigorous proportionality between oc- 
cupation number and charge density, one should, nev- 
ertheless, expect correlations between the two as far as 
general features such as oscillations or pronounced en- 
hancements/reductions are concerned. 

In view of such general correlations it may seem unclear 
why in in the cited references the charge density at the 
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FIG. 11: Oscillations in the occupation number (m) for suc- 
cessive lattice sites in the [100] direction at a band occupation 
of n — 0.8 for different coupling constants J. (Calculation for 
s.c. lattice with same Bloch dispersion as in Fig. |l|.) 

impurity is always enhanced (even divergent) indepen- 
dent of the band occupation or kp, ep etc. It seemingly 
contradicts our above argument. However, the result of 
permanent enhancement of the charge density at the im- 
purity site is non-generic and a consequence of integra- 
tions without cut-off. Integrations were performed over 
plane waves up to arbitrarily large wave vectors. Thus 
an infinite number of bands is implied with each one- 
electron state coupled to the impurity spin with prac- 
tically the same coupling strength. There is no critical 
band occupation in this effective infinite-band model be- 
cause an arbitrarily large number of electrons can stay 
at the impurity being coupled to its spin in a favourable 
way. 

IV. CONCLUSIONS 

We have provided results for the local density of 
states (LDOS) in the antiferromagnetic and ferromag- 



netic Kondo models, which has not been calculated be- 
fore with its full energy dependence. Although calcu- 
lations were performed within a simple approximation 
scheme, results prove at least in qualitative terms trust- 
worthy, allowing consistent interpretation and qualitative 
comparison with the experiment. There are clear quanti- 
tative limitations as indicated for example by the insuf- 
ficient depth of the Kondo antiresonance. Nevertheless, 
we would expect confirmation of the general features by 
quantitatively reliable methods like NRG or QMC. 

In summary, for small to intermediate couplings there 
is an antiresonance in the LDOS of the antiferromag- 
netic Kondo model (Kondo antiresonance), which has an 
analogue in the conduction-electron LDOS of the An- 
derson model. The Kondo antiresonance can be identi- 
fied with antiresonance structures observed in STS spec- 
tra of single magnetic adatoms on the assumption that 
tunneling occured mainly into the conduction band. In 
the ferromagnetic Kondo model a peak structure appears 
close to the Fermi level. It should be observable in ap- 
propriate STS experiments on adatoms with ferromag- 
netic intraatomic exchange (possibly Gd adatoms). For 
strong couplings the LDOS in both models is dominated 
by quasiparticle <5-peaks, which can be clearly related to 
the quasiparticle peaks of the exactly solvable case of an 
empty conduction band and the atomic limit. Finally, 
we have given arguments for the impurity's attractive or 
repulsive effect on conduction electrons at the impurity 
site dependent on the band occupation and for seemingly 
contradictory results gained by integrating without band 
cut-off. 
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